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NOMENCLATURE 


Z: 
x,y 
L : 
P: 
E: 
a : 


h: 

D: 


D' : 

S: 

f: 



x = (x,y) 


out-of-plane amplitude 

cartesian coordinates of plane 

typical length of vibration 

planar mass density 

Young's modulus 

Poisson's ratio 

Plate thickness 

Flexural rigidity 

D/h 

Airy stress function 

typical frequency of vibration 

x and y component of displacement from equilibrium 

Laplaci an operator 

coefficient of nonlinear term 

position vector 


i i 



x + 


- (^X» ky) : 
to a Dk 4 /p: 

(f) : 


wave vector 

angular frequency of vibration 
phase shift 


p, q, r, s: 
C: 
n : 


a : 


wave vectors 
nondimensional energy 
normalization factor 

ratios of x-side length to y-side length of plate 


f(t), P(£): 

x : 

A-j 3, C i 

Vi: 

Ri: 


fo rces 

nondimensional time 

parameter of two-mode nonlinear system 

control forces of two-mode nonlinear system 

real valued amplitudes of two-mode nonlinear system 


i i i 



CHAPTER 1 
INTRODUCTION 


Linear models have had, and continue to have, great importance in the 
simulation and control of real processes. There are, however, some systems which 
are inherently nonlinear and cannot be understood fully through linear analysis 
alone. These nonlinear systems may be encountered, for example, in fluid 
mechanics, plasma physics, and structural mechanics. 

Here, the need for undertaking nonlinear analysis arises from certain 
problems in structural mechanics. In the future, NASA expects to place into orbit 
large, extended and flexible structures; these structures may be expected to have 
oscillations large in comparison to some corresponding dimension. In this 
situation, system linearization through an assumption of only small oscillations 
is no longer possible and intrinsic nonlinearities must be retained. 

A model system for studying nonlinear structural effects is found in the so 
called von Karman equations of nonlinear elasticity. These equations describe the 
dynamics of a flat plate undergoing relatively large vibrations. It is these 
equations which will be studied here through numerical simulation by means of a 
Fourier spectral transform method. 

In this report the von Karman equations will be discussed as well as the 
numerical method for their solution. Then, preliminary results will be 
presented, both for a general and a reduced system. Finally, concluding remarks 
and possible future research directions will be given in the summary. 



CHAPTER 2 

EQUATIONS OF MOTIONS FOR A VIBRATING PLATE 


The equations of motion for a vibrating plate follow from the equilibrium 

equations for such a plate and d'Alembert's Principle. The equilibrium equations 

for a thin plate undergoing large oscillations have been determined by A. Foppl 

in 1907 and T. von Karman in 1910 and are presented in standard texts on the 
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theory of elasticity. 9 * . A few preliminary remarks are needed before 
discussing these equations. 

First, the plate will have a thickness h and the coordinate system used to 
describe its motion will be cartesian; when the plate is completely at rest, its 
neutral surface coincides with the x-y plane of the coordinate system. The 
out-of-plane vibrations of the plate will be described by the dependent variable 
1-1 (x,y,t). The following figure will be useful in discussing the equations of 
motion: 


z 



In this figure Z is a typical amplitude and L is a typical length associated with 
Z; also, f will denote a typical frequency associated with the vibration. 

The equations of motion are 
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p : mass/unit area 


E: Young's modulus, a : Poisson's ratio 

S: Airy stress function : The stress tensor is related to second spatial 

derivatives of S. 

Equations (1) and (2) are, again, derived from the Foppl-von Karman equation and 
d'Alembert's Principle; they form a set of coupled nonlinear equations which are 
accurate to a certain level of approximation. What this level is will be 
discussed presently. 
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In order to determine about how large the various terms are in the equations 
of motion, the following substitutions can be made in the various terms: 


7 ? 7 

9 9 6 - i" 2 . L 

«* p* 9x9y * • • 9 t ’ 
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Using these assignments, equation (2) gives the following order of magnitude 
estimate: S E Z 2 . The various terms in (1) can now be given order of 

magnitude estimates: 
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The first term, Tj, is not of immediate concern since it is driven by the 
other two; the ratio of the third to the second term is of interest, however, 
since it tells us the reasonableness of linearizing equation (1). This ratio is 
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T3/T2 ~ Z 2 /h 2 ; the equation may be reasonably linearized (the third term 
neglected) only if Z 2 <<h 2 . If T2 is neglected when Z 2 /h 2 is not small, 
the frequencies and mode shapes associated with a linear solution may no longer 
provide an accurate dynamical description. 


It was remarked earlier that although the equations of motion given by (1) 
and (2) were nonlinear, they were still only accurate to a certain level of 
approximation; this will be explained here. The equations of equilibrium from 
which (1) and (2) were derived can be obtained by merely replacing Tj by -P/h 
where P is the 'external' force applied to the plate. d'Alembert's Principle 
allows us to reach the equations of motion by replacing ~P by ph Z; thus, T]_ 
is really of the form 


T d 2 z a 2 z A o a 2 Z , . . a 2 z 
t i ~ p ~73 " p -.2 2 u a ax at U a u e ax ax. 

dt at a a g 


In (4) there is an implied summation over the repeated greek indices e.g., 


“1 s!/“ 2 db ■ \ 97 + % 57 Where U « a r <> the displacements in 

the x and y directions. In order to get order of magnitude estimates for the 
terms in (4) it is necessary to get such an estimate for u a . To do this, 
consider the following figure: 
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In figure 2 it is indicated that the geometry of a small section of 


vibrating plate can be replaced by a right triangle. This leads to an estimate 
of u, for Z somewhat less than L/2 » of u " Z 2 /L. Using this, along with the 
aforementioned assignments V 2 “ L -2 , etc. gives estimates of sizes of the 
various terms in (4): 
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The ratio of the second to the first and the third to the first terms are 



Although they will not be explicitly given here, there are also other terms 
of these magnitudes in (1) and (2). Thus, the equations of motion as given by 
(1) and (2) are valid as long as Z 2 « L 2 , but are not restricted by the 
particular magnitude of Z 2 /h 2 „ (If we were concerned about extremely large 
vibrations, i.e., Z 2 ~ L 2 , we would have to include these neglected terms.) 

In the study of control systems for large space structures at LaRC, the need 
to study systems for which Z 2 > h 2 has been recognized. Thus, equations (1) 
and (2) should be solved as they appear here, rather than merely solving the 
linearized versions alone. It would be informative to do this and then compare 
the nonlinear and linear solutions, both with themselves and with laboratory 
experiments. 

There are methods for solving these nonlinear equations currently available. 
There are the so-called spectral methods^ (which have been used for fluid 
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dynamical calculations over the last decade and are reaching a state of maturity) 
as well as finite difference, finite element, and Galerkin methods 5 . 

Finally, it is important to mention that in the one-dimensional case (e.g., 
the amplitude Z is a function of x, but not of y) the nonlinear term in (1) 
disappears. Remember that terms of relative size Z^/L^ have been left out of 
(1); the one-dimensional case does, in an exact formulation, still contain these 
small nonlinear terms. For our purposes, however, and to a high degree of 
accuracy, we may view the one-dimensional case as inherently linear while the 
two-dimensional case will be, in general, nonlinear. 
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CHAPTER 3 

SPECTRAL METHOD OF SOLUTION 


The equations of motion for a vibrating plate (the von Karman equations) 
are, again. 


2 2 2 2 2 
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O' = E : Young's modulus, a : Poisson's ratio 

12 ( 1-a 

h: plate thickness 

S: Airy stress function : The stress tensor is related to second spatial 

derivatives of S. 

f: external force applied to plate 


The applicability of these equations has been discussed in the previous chapter. 
The basic results were: (a) that the out of plane amplitude could be 

considerably larger than the plate thickness (in which case the nonlinear terms 
in (5) become important) and (b) that the amplitude of a particular oscillatory 
solution was considerably less than the wavelength associated with that 
solution. 


If the nonlinear and external force terms are neglected in (5), we have a 
linear equation of motion (which is sufficient in itself only for general 
one-dimensional motions or very small amplitude two-dimensional motions.) The 
linear equation has a general solution of the form 
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(7) 
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where x = (x,y), £ = (k , k ), = D k^/p and the constants A, B, C» D, and <j> 

x y 

are determined by the boundary and initial conditions. If the non-linear terms 
are included in (5), then an analytic solution such as (7) is, in general, not 
obtainable. 


In the nonlinear case, an alternative method of solution is to proceed 
numerically. There are several numerical methods to choose from, for example, 
finite element, finite difference or spectral (Galerkin). Here, we will discuss 
the spectral method, in particular, on the (Fourier) spectral method of Orszag 

ft 1 

and Patterson * . 


In this spectral method the dependent quantities in the equations of motion, 
(5) and (6), are expressed in terms of Fourier series: 
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where Z (£, t) 

= Z * 

'(-£, t) and S (£, t) = 

S *(-£, t) (*: 

compl ex 

conjugation) , 

N and M 

are the number of points 

in the x and y 

directions , 

respectively; 

the summation limits will be discussed in Chapter 4. Whether 


we are discussing physical variables Z(x), S(x) or Fourier modes Z(t)» S(t) is 
determined by their respective arguments (often the time argument will be 
omitted.) This choice implies that we are imposing periodic boundary conditions; 

8 

B.C.'s can be discontinuous and Lanczos* Tau method can handle arbitrary B.C.'s . 
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Notice that these boundary conditions keep the trigonometric spatial 
eigenfunctions of the linear solution (7) while eliminating the hyperbolic ones, 
also iri (7), from the nonlinear solutions (8) and (9). This is not a fundamental 
limitation since the trigonometric functions still provide a complete orthogonal 
basis. 

The time dependence of Z (x, t) and S (x, t) is now contained in the time 

dependence of the Fourier coefficients Z (t, t) and S (t, t). The equations 
which describe the time evolution of the Fourier coefficients are determined by 
substituting (8) and (9) into (5) and (6); after some algebraic manipulation we 
have 

P d 2 z(t, t) + D . k 4 z( t, t) - 4= l |5 *3| 2 Z (P. t) S (q, t) = f (t, t) 

dt /NM X ' ' (10) 

p+q 

S (q,t) E |r x si 2 Z(r, t) Z(s, t) (11) 

2q 4 /M Jq 1 1 

r+s 

Here f($, t) is the £ - Fourier mode coefficient of the external force and 
all t = 0 coefficients are omitted from the solution. (The components of all the 
modal vectors 1c, p, q, r, and s appearing in (8), (9), (10), and (11) are 
integers, while the components of x in (8) and (9) take values between 0 and 2n.) 

In the Orszag-Patterson spectral method, an appropriate time-differencing 
scheme is chosen to advance Z($, t) in (10) forward in time. The essence of the 
method is not to solve for the nonlinear terms in (10) and (11) by simply 
summing them in t-space (which takes a very long time), but instead the method 
finds the coefficients of the derivatives of S(x) and Z(x), uses an FFT (Fast 
Fourier Transform) to transform to x - space, forms the nonlinear terms of (!5) and 
(6) there, then on a discrete grid of points, say N x M, where N and M are powers 
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of two. The x and y components of x take on the values x n = 2tt n/N and y = 

2 it m/M 0 £ n £ N-l and 0 £ m £ M-l). The ratio in speed between the transform 

method and the direct sum evaluation is (logN)/N where N is also the number of 

distinct k or k values (i.e., the number of grid points in x - space and t - 
x y 

space are equal). Another important part of the method is to properly de-alias 

the nonlinear coefficients determined by the transform technique; this is 

6 7 

discussed by Orszag , and Patterson and Orszag . 

The purpose of this chapter has been to briefly describe the Orszag-Patterson 
spectral method as applied to the problem of nonlinear vibrations of a thin 
plate. The method can be used to study the importance of nonlinearities, 
particularly as the mean amplitude to plate thickness ratio is increased. This 
could demonstrate the need to more closely consider nonlinear effects in our 
attempts to understand and control the dynamics of large space structures. In the 
following chapters the implementation of this method and some preliminary results 
will be discussed. 
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CHAPTER 4 


SPECTRAL EQUATIONS: FURTHER CONSIDERATIONS 

The spectral equations developed in the last chapter can be put into a 
nondimensional form suitable for studying the nonlinear vibrations of a 
rectangular plate. The end result of doing this is shown in Figure 3; it will 
now be assumed that the number of points in the x and y directions will be the 
same and equal to N. 

The equations in Figure 3 can be consolidated to eliminate the Airy stress 
function (i.e., its Fourier coefficients). The result of this action is 
presented in Figure 4a; the coefficient of the cubic nonlinear term has been 
explicitly symmetrized. The reason for this is the following: First, the time 

dependent equation in Figure 4a is multiplied by Z*(£, t) and the resulting 
equation is added to its complete conjugate; then, using the fact that N (£, p, 
r, s) is completely symmetric in its arguments allows the time dependent 
equations in Figure 4b to be arrived at. 

As is evident in the equation for d£/dx in Figure 4b, in the absence of an 
applied force, the quantity 5 is constant (but only for an infinite Fourier 
expansion). This conserved quantity E, is the energy of the vibrating plate; the 
term with coefficient 0 is the 'nonlinear' contribution; if 0=0, then ? = 
is the 'linear energy'. There does not appear to be any other constant of the 
motion, as is the case in analogous representations of fluid and magneto-fluid 

9 

systems . 

The behavior of the energy for a system represented by a truncated (finite) 
Fourier representation will be examined in the numerical example of the next 
chapter. Let us conclude the present chapter with some remarks concerning the 
contribution of the t = 0 modes (coefficients) and on the statistical nature of 
the nonlinear dynamics. 

Careful consideration of the form of N (£, p, r, s) shows that as t -> 0, 

N (...)>0, regardless of the values of p, r, and s„ Thus, the \ = 0 Fourier 
coefficient has no effect on the dynamics of the plate and corresponds only to 
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the overall velocity of the chosen inertial reference frame. In the numerical 
work that follows, the zero mode coefficient will always be set equal to zero. 

As in the case of nonlinear fluid systems, the presence of nonlinear terms 
introduces random and perhaps chaotic behavior into the time evolution of the 
Fourier coefficients. This randomness may be expected to manifest itself as a 
short correlation and autocorrelation time of the Fourier coefficients; this, in 
turn, allows for a statistical description of the dynamics (el astodynamic 
turbulence). This description can be effected through the use of a canonical 
(Gibbs) distribution* 0 where the distribution function is 

D U) ■ IJULhzil (12) 

Here a is a parameter (inverse temperature) which is chosen so that D (?) is an 
extremum for £ = l , the actual energy of the system. The quantity n is a 
normalizing factor: 

n = / exp (-<*£■) | d l R (t) d V R (t) d Zj(t) d Vj(tc) (13) 

Is 

where Z R j(t:) and V R j(£) = d Z R j(£)/dx are the independent coordinates of the 
statistical phase space (Z R j(t): real, imaginary parts of Z(t)). 

In fluid dynamical systems a similar statistical description can be utilized 

because constants of the motion remain constant even when the Fourier series are 

g 

finite, as in the case in any numerical simulation. Such is not the situation 
for any finite expansion of an el astodynamic system, however, although certain 
model systems can be developed in which the energy remains constant even with a 
finite number of modes. This topic will be touched upon again in the following 
chapters. 
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Z = RMS vibrational 
amplitude 



E: Young's Modulus a: Poisson's ratio D' = 

p: density h: plate thickness 

P(j<): Fourier coefficients of nondimensional force 
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Fig. 4a Consolidated Equations of Motion 
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Fig. 4b. Constant of the Motion: the Energy 
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CHAPTER 5 

NUMERICAL SIMULATION 


The numerical solution of the nonlinear von Karman equations has attracted a 
good deal of effort since their inception, as is evidenced by a survey of the 
open literature. Although spectral (Galerkin) methods seem to be widely applied, 
the efforts usually end in a linearization of the problem. An example is found 
in the article by Crawford and Atluri^ entitled "Non-linear Vibration of a Flat 
Plate with Initial Stress". Crawford and Atluri present a solution by means of 
Galerkin approximation of the von Karman equations using the (linear) 
eigenfunctions appropriate for simply supported boundary conditions 
(two-dimensional sine series). They specialize to a problem containing only four 
modes (i.e., four independent sine series coefficients) and use a perturbation 
technique called "multiple time scaling" to solve the nonlinear ordinary 
differential equations for the sine series coefficients. They then show that 
their results are commensurate with prior perturbative analysis. 

The intent here is to get away from perturbative analysis, wherein one 
method may easily be consistent with another, and present an alternative 
approach. This new approach is the Fourier spectral transform method discussed 
in Chapter Three. In this chapter, the numerical implementation of that method 
will be discussed. In particular, pertinent computational details will be 
considered, some preliminary results presented and practical limits on computer 
resource utilization will be given. 

In order to create a computer code to simulate nonlinear dynamics using a 
Fourier spectral transform method, some essential features must be incorporated. 
First, there has to be a computer system available suited to intensive vector 
computation. At the NASA Langley Research Center, the CDC Cyber 203 (STAR) is 
such a machine (although it does have its inherent 1 imitations--to be considered 
shortly). Next, there must be a reasonably quick and efficient Fast Fourier 
Transform (FFT) subroutine resident on the system (the Cyber 203 mathematics 
library contains the appropriate routine: Q4FFT2DR). The FFT subroutine is 

important because a majority of the computational time is spent performing 
transforms. 
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In the computer code, two important considerations are the time integration 

scheme and the de-aliasing procedure. In the preliminary code developed here, a 

12 

second order 'two step' method was utilized; this second order method was used 
rather than a higher order one so as to reduce the number of arrays needed. A 
de-aliasing procedure is needed because the nonlinear interaction, by its nature, 
feeds energy into both higher and lower order modes from those modes which have 
been initially excited. Thus, when FFT’s are performed at each time step, it is 
possible to get aliased values for those modal coefficient which are being 
retained. This problem is circumvented by either performing 2N point transforms 
and keeping only the first N (or less) coefficients, or by using the 
advanced methods of Patterson and Orszag 6 ’^; in this preliminary development the 
first procedure was used. Specifically, only those modes Z (t) w|t| = k _< K were 
retained; all modes outside the circular region defined by k = K were reset to 
zero after each forward FFT. 

As has been mentioned, the numerical technique used here has been 

g 

effectively utilized to study problems in nonlinear fluid dynamics . In these 

fluid applications the modal dynamic equations contain quadratic non-linearities, 

while in the el asto-dynamic case the nonlinearities are cubic. Another essential 

difference is that the fluid equations contain at most second order spatial 

derivatives while in the von Karman equations there is the biharmonic term, which 

is of fourth order. Therefore, if K is the largest wave vector magnitude in both 

the fluid and el astodynamic simulations, then the time step 
-2 -4 

sizes scale as K and K , respectively, in these simulations; thus, while 
numerical simulations of fluid turbulence (nonlinear dynamics) to a certain 
(moderate) level of resolution may currently be feasible, the same level cannot 
be reached in elasto-turbulence simulations (unless a huge computer resource 
investment is made). 

In the preliminary testing and evaluation of the transform code for the 
nonlinear von Karman equations the plate was represented by a grid of 16 x 16 
points. The k-space cutoff was K=4; the dimensionless time step size was 
(2048)”*, as any lower values allowed instabilities to grow. The initial 
conditions established for the Fourier coefficients F(t,t) were |F(t,0)| 
exp [(jlc L ) /2 s^ 3 where fc = 2 and the standard derivation was s = 2. The 
phases of the coefficients were chosen randomly and the overall energy was 
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normalized so that the 'linear 1 part was one (i.e., that part of the energy ? 
in Fig. 4b which remains when $ = 0). 

With these initial conditions the code was run with 3=0 (linear) and 3 = 
256 (nonlinear). For 3 = 0 the energy in each mode remained constant during the 
time evolution of the motion (as expected.) For 3 = 256 modal energy did not 
remain constant and some modes were preferentially enhanced; this is illustrated 
in the sequence of Figure 5. Figure 5a is the linear part of the modal energies 

? L (£) at t = 0 (for 3=0 this energy distribution did not change with time). 
Figures 5 b and c show the transferral of energy for 3 = 256: energy in the mode 

£ = (3,0) and (-3,0) is augmented by nonlinear processes. 

For a 1 cm thick square aluminum plate of side length 3 meters, the time of 
evolution shown in Figure 5 is 0.86 seconds; for a side length of 30 meters, this 
time becomes 86 seconds. The frequencies corresponding to these two cases are, 
for a side length of 3 meters, lowest (|£| = 1): 11.56 Hz and augmented (£ = 
(3,0)): 105 Hz; for a side length of 30 meters, these become, lowest: 0.116 Hz 
and augmented: 1.05 Hz. There were no specific boundary conditions imposed; 

this was a 'homogeneous' el asto-turbul ent simulation. 

4 

The presence of terms of order K in modal von Karman equations 
(Figure 3) requires the investment of relatively large amounts of computer time. 
Table 1 presents some pertinent parameters for a typical low resolution run. 

Table 1. Computational Parameters 

Grid Size: 16 x 16 

k-space cutoff: K = 4 

Total time: x = 100 

STAR cpu-sec/time step: 0.048 sec 

Total cpu time: 2.78 hours 

Number of retained modes: 50 

As can be seen from Table 1, a good bit of time must be spent to perform a 
nonlinear el astodynamic simulation. So much time that the resources available 
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for this purpose were quickly expended, while interesting parts of the simulation 
remained untouched. In Figure 6 is a comparison of the behavior of the total 
energy £ and the linear energy between times x = 0 to 120. Here, it is 
evident that the total energy is not constant, as would be the case if an f 
infinite number of modes could have been utilizied; instead, the energy falls by 
about 1 % of its initial value. This 1% can be thought of as 'heat* loss: energy 

has been cascaded over the = K cutoff and is no longer dynamically useful. 
Note that the nonlinear contribution to the energy, 5^ = £ - $ L , also feeds the 

linear energy and causes it to rise temporarily (although the total energy 5 
monotonically decreases). 

The code described herein is a fully nonlinear implementation of the von 
Karman equations, subject only to discretization error. If the computer 
resources are available, there are many applications which can be profitably 
studied, including flutter analysis, large space structure (LSS) dynamics and 
non-linear control; appropriate boundary conditions can be applied using Lanczos' 
Tau method. However, because the necessary resources are not, at present, 
available, a multi-mode, high resolution simulation is not feasible. This forces 
us to use a low-resolution model in which only a few modes are present. Such a 
model will be studied in the next chapter: the nonlinear dynamics of a von 

Karman system restricted to two modes will be investigated. 
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CHAPTER 6 


DYNAMICS AND CONTROL OF A MINIMAL 
NONLINEAR SYSTEM 

Dynamics 


The results of the previous chapter have shown that if a truncated Fourier 
representation is used to describe the amplitude Z (x, t) , and if the retained 
Fourier modes are initially arbitrarily excited, then the energy of the model 
system is not necessarily conserved. These results also illustrated the high 
computational cost incurred when retaining even a moderate number of Fourier 
coefficients. Since the purpose of this report is to gain a perspective on the 
control of a nonlinear elastic system, it would be useful to have a model 
numerical system which was inexpensive to run and in which numerical dissipation 
did not cloud the effects of the control scheme under study. Just such a model 
system will be discussed here. 

Consider the modal equations represented in Figure 4a; suppose that there 


are only two independent Fourier modes initially excited: Z^ = Z(t^ ,t) and Z^ = 

Z (t 2 ,T) and that all other modes begin and remain unexcited. The model system 
consistent with these assumptions has the form 


2 ^ 

r- + < A 1 + c | z 2 | z > 



(14) 


d 2z 2 


2 + ^2 + C Z^ ) Zj = Ug 


dx 


From Figure 4a we have the correspondence 


Aj = </ (&j, a) 


C = 


i |^1 X + ^2 1 4 + |^1 “ ^ 


-4 


1-4, 


3N 
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In order to study the system of equations (14) the explicit values of e, N» tj, 

"L and a are not needed, so long as the A- and C are non-negative. It should be 
^ J 

noted that (14) are in general complex; however, since Z (t, t) = Z*(-t,t), 
the complex conjugates of (14) are not independent of (14) themselves, and the 
two (complex) equations of (14) are thus the complete set of system equations. 

Since the Z i are generally complex, the two explicit equations in (14) are 
implicitly four. However, if the Z^ are assumed to be each either fully real 
or fully imaginary, then there are only two equations; this then is the smallest 
reduction of the von Karman system of equations which remains nonlinear: two 

equation in two unknowns. 

These two, coupled nonlinear equations can represent a real physical system 

if Zj = i Rj, Z 2 = i R 2 (i = /^T ; R ^ : real). The system then corresponds to a 

simply supported vibrating plate, the situation considered by Crawford and 
Atluri. The system of equations (14) becomes 

d 2 Rl 2 

— y~ + (A x + C R 2 ) Rj = Uj (15a) 

d 2f L 2 

— Y- + (A 2 + c RJ) R 2 = U 2 (15b) 

&x L 

This basic system can now be studied to gain some insight into the nature of a 
more complex multi-mode system. 

To begin, if (15a) is multiplied by dR^dt and (15b) by dR 2 /dT, there 
results the following two equations: 
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(16a) 


j | dR* a i o 1 o ^ ^ 1 dR- 

clT C 7 ( Tr } + 7 A 1 R l ] + 7 c r 2 H7" = U 1 H7~ 


2 

i i (IRa 9 1 9 19 ^^9 ^Rq 

dT ^7 ( dT" } + 7 A 2 y + 7 c R l dT" = u 2 dT" 


(16b) 


If this system is linear, C = 0; furthermore, if there is no forcing, 

U 2 = 0 and the following quantities are individually conserved (k = 1,2) 


1 2 2 

« k - z r - ( dT> + \ v 


(17) 


The £ k (k = 1,2) are the 'linear' energies of the two modes. Now, suppose 
C * 0 but U 1 = U 2 = 0; and £ 2 are no longer individually conserved, but a 
single conserved quantity still exists, for if (16a) and (16b) are added together 
the result is 


d 

dx 



+ «2 + 2 - C ^ 


R*> 


dR, 
J 1 dT 


dR. 


+ U. 


2 dx 


(18) 


If the modal control forces and U 2 are zero the total energy of the system £ = 
12 2 

= ?1 + i; 2 + 7T CR^ R 2 is obviously conserved. 

The system of equations (15a, b) can be easily solved numerically, this was 
done using the same time-stepping scheme as for the multimode system 
(Ax = 0.01) although here the nonlinear terms are evaluated directed (no 
transfrom). This was done using the same time integration scheme as for the 
multi-mode case, with a dimensionless time-step of 0.01; here, the nonlinear 
terms were evaluated directly, rather than by transform. As an example, let 
= 1.0, l \2 -■ 2.0 and C e {0, 1/4, 1/2, 1, 2, 4, 8, 16}; this procedure 
allows the behavior of the system to be studied as the nonlinearity is 
increased. The results are shown in Figures 7a-h; notice that although 5 is 
conserved, and ^ are not, as 6 increases in value. Energy is apparently 
being transferred between the two modes through their mutual non-linear 
coup! ing. 
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In Figure 8 another interesting aspect of the dynamics is revealed; there, a 
comparison is given between the temporal behavior of the total energy £ and the 
sum of individual linear energies and Notice that ^ * £ in general; 
i.e., there also exists a 'non-linear* energy of coupling between the modes which 
fluctuates between zero and roughly half of the total energy. In Figure 8 we see 
that the non-linear coupling energy goes through many fluctuation cycles as 
energy is being shifted back and forth between modes. In contrast consider 
Figure 6; there the multi-mode system has apparently not even gone through one 
fluctuation cycle, suggesting that most of the interesting dynamical history 
remains to be explored (given sufficient computer resources). 

Control 

The minimal nonlinear system described in this chapter was used to study the 

effect of a linear control law on a structure which was inherently nonlinear, as 

this situation may occur in practice. The discrete control law chosen was that 

13 

of Williams and Montgomery which uses the two most recent past values of 
displacement and control inputs to calculate the new control inputs at each time 
step. The fact that this controller uses two past values of displacement implies 
the presence of rate feedback; this is important because it is evident from (18) 
that modal control inputs opposite in sign and proportional to the instantaneous 
modal rates will always decrease the total energy (the 'Lyapunov' function of the 
system) . 

The characteristic equation of the system with the controller of Williams 
and Montgomery has four roots. Here, it was assumed that these were two pairs of 
conjugate roots; the real and imaginary parts of one pair were both equal to 
0.3535 and +_ 0.3535, respectively, while the other pair had real and imaginary 
parts of 0.4949 and 0.4949, respectively. This placed the roots well inside 
the unit circle in the complex plane, ensuring stability. Using this control 
scheme, the nonlinearity parameter C in (15a, b) was set to a number of values, 
starting with zero, to examine the effect of increasing the nonlinearness of the 
system on the controller. The sampling time was, in dimensionless units, 0.20. 
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The behavior of the amplitudes and R£ of (15a, b) are shown in Figures 9, 
10, 11, and 12. In Figure 9 and 10, the system parameters are = 1, = 2.25, 

C = 0: the system is linear, the modes are uncoupled and damp out smoothly. In 

Figures 11 and 12, = 1, A 2 = 2.25 and C = 8: the system is highly non-linear, 

the behavior is initially erratic, but still damps out, though not as quickly as 
when C = 0. Notice that during the initial erratic phase, the oscillations of 
both R.| and R ? do not have zero mean. (Figures 9 and 10, the ideal responses 
of a linear system, should be compared with Figures 11 and 12, respectively, the 
actual responses of the nonlinear system.) 

This result indicates that although the dynamic response of the nonlinear 
system to linear control is erratic at first, the system is ultimately 
controlled. It would be interesting to study this phenomenon further, 
particularly on the multi-mode system discussed in the previous chapter. In so 
doing, the presence of spillover could be investigated, a problem which is 
explicitly avoided in the minimal nonlinear system. The effect of nonlinearities 
on desired controller response, along with other interesting topics, could also 
be examined (provided the computational resources are available). 
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CHAPTER 7 
SUMMARY 


In this report the nonlinear dynamics and control of a vibrating rectangular 
plate were investigated. The dynamics of the plate were simulated by solving the 
von Kannan equations through the use of an Orszag-Patterson type Fourier spectral 
transform method. Although arbitrary boundary conditions could be implemented by 
using Lanczos' Tau Method, this was not done due to the lack of computer 
resources. 

In order to overcome the computational cost inherent in a full-blown, 
multi-mode simulation, a reduced order model was studied. This reduced model was 
a minimal von Karman system: two coupled nonlinear differential equations in two 

unknowns. This minimal set contained all the qualitative details of a multi-mode 
simulation: nonlinear coupling and the presence of a Lyapunov function. 

Preliminary results indicate that controllers with at least implicit rate 
feedback, even if predicated on a linear modeal , can eventually control and 
dampen the nonlinear oscillations, as long as the nonlinear terms are not too 
dominant. 

The minimal and the multi-mode systems discussed herein can be profitably 
used to study a number of interesting problems. Different control laws can be 
studied using either the minimal or multi-mode system as the plant; in the 
minimal system, spillover is explicitly avoided, while in the multi-mode system 
its full effects can be investigated. Lyapunov theory can also be applied in 
more detail to these nonlinear systems to help shed light on questions of 
control . 

Additional problems of interest include the relation between modal 
autocorrelation time and control system sampling rate, and the phase space 
structure of both the conserved and damped nonlinear dynamic system. (Other 
nonlinear systems have already been shown to possess interesting and unexpected 

1 c 

dynamic behavior ). 
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The end to which this research is directed is, of course, to understanding 
and controlling the dynamics of real-world systems. Amongst these are large 
space structures, fluttering wings and other extended, relatively flimsy 
structures of interest to NASA. As these structures are difficult to 
experimentally simulate on the earth's surface, and since we would like to know 
as much as possible before putting these objects into orbit, a numerical 
simulation provides a viable alternative with which to gain the requisite 
knowl edge. 
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